Information propagation for interacting particle systems 
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We show that excitations of interacting quantum particles in lattice models always propagate with a finite 
speed of sound. Our argument is simple yet general and shows that by focusing on the physically relevant 
observables one can generally expect a bounded speed of information propagation. The argument applies equally 
to quantum spins, bosons such as in the Bose-Hubbard model, fermions, anyons, and general mixtures thereof, 
on arbitrary lattices of any dimension. It also pertains to dissipative dynamics on the lattice, and generalizes to 
the continuum for quantum fields. Our result can be seen as a meaningful analogue of the Lieb-Robinson bound 
for strongly correlated models. 



How fast can information propagate through a system of 
interacting particles? The obvious answer seems: No faster 
than the speed of light. While certainly correct, this is not 
the answer one is usually looking for. For instance, in a clas- 
sical solid, liquid, or gas, perturbations rather propagate at 
the speed of sound, which is determined by the way the par- 
ticles in the system locally interact with each other, without 
any reference to relativistic effects. We would like to under- 
stand whether a similar "speed of sound" exists for interact- 
ing quantum systems, limiting the propagation speed of local- 
ized excitations, i.e., (quasi-)particles. For interacting quan- 
tum spin systems, such a maximal velocity, known as the 
Lieb-Robinson bound Offil, has indeed been shown. While 
it seems appealing that there should always be such a bound, 
systems of interacting bosons can show counterintuitive ef- 
fects, in particular since the interpretation of excitations in 
terms of particles is no longer fully justified; in fact, an exam- 
ple of a non-relativistic system where bosons condense into a 
dynamical state which steadily accelerates has recently been 
constructed [5|. This example suggests the disturbing possi- 
bility that our intuition is wrong, and only relativistic quantum 
theory can provide a proper speed limit. 

There are many important reasons, both theoretical and ex- 
perimental, to investigate information propagation bounds in 
interacting particle systems. It turns out that such bounds lead 
directly to important, general results concerning the cluster- 
ing of correlations in equilibrium states J2). Lieb-Robinson 
bounds facilitate the simulatability of strongly interacting 
quantum systems — the mere existence of a Lieb-Robinson 
bound for a quantum system can be used to develop general, 
efficient, numerical procedures to simulate the dynamics of 
lattice models (SJ. From a more practical perspective, new ex- 
periments allow one to explore the non-equilibrium dynamics 
of ultracold strongly correlated quantum particles — bosonic, 
fermionic, or mixtures thereof — in optical lattices with un- 
precedented control [7 8]. In such experiments, it is impor- 
tant to understand how the particles move: For example, when 
studying instances of anomalous expansion, it is far from clear 
a priori whether it is possible to identify a meaningful speed 
of sound at all. 




FIG. 1. Schematic representation of the "light cone" of particles ini- 
tially placed into a region R of a lattice (yellow circles) and then 
propagating in time t in a way governed by an interacting quantum 
model, outside of which the influence of these particles is exponen- 
tially suppressed. 

The original Lieb-Robinson bound already applies in a very 
general setting, namely, to any low-dimensional quantum spin 
system, and to any fermionic system confined to a lattice. It 
is therefore tempting to extend the original argument to other 
settings, in particular, to systems of interacting bosons; unfor- 
tunately, all attempts to do so have run into insuperable diffi- 
culties for systems with nonlinear interactions, including the 
Bose-Hubbard model. The reason for the failure of the orig- 
inal Lieb-Robinson argument is fundamentally connected to 
the unboundedness of the creation operator for bosons: The 
Lieb-Robinson velocity depends on the norm of the interac- 
tion, which is unbounded for, e.g., bosons hopping on a lattice, 
and examples without a speed limit can be constructed |5 1. 

In this Letter, we show how these difficulties can be over- 
come by considering the right question concerning the prop- 
agation of information. Our approach allows us to determine 
Lieb-Robinson type bounds for the maximal speed at which 
information can propagate through systems of interacting par- 
ticles in a very general scenario: In particular, it applies to 
systems of interacting bosons, as well as to fermions, spins, 
anyons, or mixtures thereof, both on lattices and in the con- 
tinuum. Moreover, it can also be applied beyond Hamiltonian 
evolution, such as to systems evolving under some local dis- 
sipative dynamics. 
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The type of system we have in mind is exemplified by the 
Bose-Hubbard model, a model of bosons hopping on an arbi- 
trary lattice G of any finite dimension and interacting via an 
on-site repulsion, 



h.c. 
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where the first summation is over neighboring sites on the 
lattice, bj is the boson annihilation operator for site j, and 
hj = hjbj is the number operator. The natural distance in 
the lattice will be denoted by d(-, •), e.g., d(j, k) = \j — k\ 
for a one-dimensional chain. While we will, for clarity, focus 
our discussion on the Bose-Hubbard model, our arguments di- 
rectly generalize to models of the form 
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where the b St j are annihilation operators for bosons, fermions, 
or even anyons of species s = 1, . . . , S at site j, and n s j — 
P s j b s ■ ; the species could for instance refer to an internal spin 
degree of freedom. The interaction between the particles is 
characterized by / which can be an arbitrary function of the 
local densities, and may involve higher moments of the par- 
ticle number, or even non-local interactions. Moreover, our 
argument also applies to time-dependent Hamiltonians of this 
form, as long as the tunneling amplitude r() is bounded. 

The scenario we consider is described by the Bose-Hubbard 
model on a lattice G, where in the initial state all sites are 
empty (i.e., (hj) — 0) except for the sites in a region R which 
can be in an arbitrary initial state with finite average particle 
number. Note that the region R may very well encompass the 
major part of the lattice. What we are interested in is how 
fast these bosons will travel into the empty part G\R of the 
lattice, as a function of the distance d(-, •) on the underlying 
graph. In particular, we would like to find a "speed of sound" 
for the bosons, that is, a velocity v such that for any region S 
in G\R with d(S, R) > I [i.e.: d(s, r) > I Vs G S, r e R], 
and for all times t for which vt < I, the expectation value 
of any observable Os on S is equal to the expectation value 
of the vacuum, up to a correction which decays exponentially 
away from the light cone, e 7 ^* - ^. 

To start, we consider the Bose-Hubbard model Hbh an d 
focus on measurements of the local particle number operators 
hj. This corresponds to looking for bosons at the initially 
empty sites, and thus captures the most natural notion of par- 
ticles propagating into a region. Let us denote the initial state 
by p(0), which evolves according to 

p(t) = -i[H BH ,p(t)} 

for t > 0. As we are interested in the speed at which particles 
in the Bose-Hubbard model propagate, let us try to understand 
how the local particle densities 

aj(t) = tr(hjp(t)), jeG, 



evolve under Hbh- To this end, we derive a bound on the rate 
at which ay (•) changes, which in turn leads to a bound on the 
velocity at which particles can propagate through the system. 
It holds that 



= -«tr(fij [H B n,p(t)]) = -itv([hj,H Bli ]p(t)) 
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where the summation runs over all sites k neighboring j, 
d(j, k) = 1. Since we are only interested in an upper bound 
on this rate of change, we now consider |ou(£)| and apply the 
triangle inequality to obtain 



\d j (t)\<2rJ2Hb{b j p(t))\ 



(4) 



To bound this term we use the operator Cauchy-Schwarz in- 
equality, viewing 

tr(b{b jP (t)) = (b k p^(t)rb jP ^(t)) 

as a Hilbert-Schmidt scalar product of & J p 1 / 2 (i) and 
6fep 1 / 2 (t), where p 1 / 2 ^) is the matrix square root of pit). 
This gives rise to 



tr(S£V(t))| < {tr(blb k p(t))tr(b]b jP (t)) 
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Combining this with Eq. (|4]), we obtain a set of coupled dif- 
ferential inequalities 



ij(t)l < 2tJ2 
U,k) 
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which, using ^fxy < (x + y) /2, yields the linearized system 

\a j (t)\<TiVa j (t) + J2o'k(t)j , 
V <i,k> J 

where T> is the maximal vertex degree of the interaction graph. 

We are interested in the worst-case growth of otj (t ) as t pro- 
gresses. This will occur when we have equality in the above 
expression (i.e., the derivative is as large as possible), and thus 
a bound 7fe(t) > ctk(t) is given by the solution of the linear 
system of differential equations 



7j(t)=r(Z>7j(t) + X)7k(*) 



which fulfills 7j(0) = Oj (0). This solution has the form 

7(*) = e^e^7(0), 

where M is the adjacency matrix of the lattice, i.e., Mj t k = 1 
if d(j, k) = 1 and otherwise, and 7 := (jk)keL- This yields 
an upper bound 



a(t) < e 



VrtrMt -> 
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for the expected particle number at time t for any site, for 

a := (a k ) keL . 

In order to understand how quickly particles propagate 
from the initially occupied region R into a region S with 
d(R, S) > I, we need to consider the off-diagonal block of 
e Vrt e TMt correS p 0nc ji n g to those two regions. Thus, in order 
to obtain a light cone with an exponential decay exp(vt — I) 
outside it, we need to understand how rapidly the off-diagonal 
elements of the banded matrix M grow under exponentiation 
e rMt . This can be done by applying Theorem 6 from Ref. J5J, 
which yields for the (i, j)-th element of exp(rMi) the bound 

[exp(TMt)] id < C e «ot-d(i,i) 

with velocity vo = xAr, where x ~ 3.59 is the solution of 
xlriX = X + LA = IIMHoo/2 depends on the lattice dimen- 
sion, and C = 2\ 2 /(x — 1) ~ 10. Together with the prefactor 
exp(2?-rf), this gives a Lieb-Robinson velocity v = vq + Vt 
ifTUl . For the scenario of an empty lattice with particles ini- 
tially placed in a region R, this implies that for any j with 
d(j,R)>l, 

Oj{t) < Ce vt ~ l «fc(0) = CN Q e vt ~ l , (6) 

k£R 

i.e., up to an exponentially small tail, the particles propagate 
with a speed no faster than v, independent of their initial state. 
Here, Nq — J^keR a k(fy — (N) is the total number of parti- 
cles in the system (i.e., the expectation value of the total par- 
ticle number operator N — J2j Note that while this (un- 
surprisingly) means that the strength of the signal observed 
may depend on the number of bosons initially put into the 
system, the maximum propagation speed v does not depend 
on Nq. In fact, for a purely harmonic one-dimensional model 
for C/ = 0, the exact speed of sound is indeed linear in r, so 
the above bound is tight up to a small constant prefactor. 

Having understood how to obtain a bound on the propa- 
gation speed of particles, we now turn to more general ob- 
servables. First, let us show how we can bound the higher 
moments of the particle number operator. For p > 1, 

a? (t) = tr(f$p(t)) 

= 1 ^2ix(h j h^- 1 p N p(t)P N ) 

N 

<J2 tT ( fl 3 NP ~ lp Np(t)PN) (7) 

N 

8 ^2 NP ^ (CNe vt ~ l ) tr(p(t)) 

N 

= C (N p ) e vt ~ l , 

where Pn projects onto the subspace with a total of N par- 
ticles, and we have used that Eq. |6} applies to each sub- 
space with fixed particle number independently as the Hamil- 
tonian commutes with Pjv- Here, (N p ) denotes the (time- 
independent) expectation value of the p-th moment of the total 



particle number operator. This proves a Lieb-Robinson bound 
for the higher moments of the particle number operator. 

Let us now turn our attention towards arbitrary local ob- 
servables Aj. Any such observable can be written as Aj = 
J2 P q c p ,q{b^) p b q j, and we have thus that 

|tr(i iP (t))|<^|c M ||tr[(StfS>(i)]| (8) 

<El^l (tr[(6t)"6jp(t)]tr[(St)96« p (t)]) 1/2 . 

In turn, for p > 0, 

tr[(6t)*fi£p(t)] = trfofo - 1) • • • (hj -p+l) P {t)} 
p 

^Vft^C/-' (9) 
r=l 

by virtue of Eq. (rob, for some constant C p . If p = 0, we 
trivially have tr[p(ij] = 1. Together, this yields abound 

\tr(A jP {t))\ < C'e vt ~ l 

if c o,q — c p.o — for all p and q, and 

|tr(i i/0 (t))| < C'e {vt - l)/2 

otherwise, where we have assumed that \ c p,q\ i s finite, and 
used that w.l.o.g. co.o = 0. In both cases, this means that 
outside the light cone given by vt — I, tr(Ajp(t)) decays ex- 
ponentially; however, the decay is on double the length scale 
in the latter case. 

Finally, observables acting on more than one site can be 
bounded analogously to the local case: Any two-site operator 
acting on sites j, k can be written as the sum of terms AjA k , 
and 

\tr{A 3 A kP {t))\ < (tv(A]A jP (t))tv(A k Alp(t))y /2 . 

The terms on the r.h.s. are local observables which can be 
bounded as before by exp(vt — I), yielding the same expo- 
nential bound for two-site — and recursively for many-site — 
observables. (Note that there exist cases where terms which 
are bounded by cxp[(vt — l)/2] only appear, and in addition 
one of the A's above could be the identity. Thus, bounds of 
the form exp((ui — I) /k) can occur, where k can grow expo- 
nentially in the block size. This, however, still implies that the 
signal is exponentially small outside the light cone.) 

While we have illustrated our arguments for the Bose- 
Hubbard model, they generalize straightforwardly to the more 
general class of models described by Eq. Q. First, it is clear 
that we can replace the on-site replusion and chemical po- 
tential in the Bose-Hubbard model by any type of interac- 
tion (even a non-local one) which only depends on the par- 
ticle numbers, since any such term vanishes in the commu- 
tator [hj , H] in Eq. (13b. Second, for systems that contain 
several types of bosons the same arguments apply: Such sys- 
tems can be modelled using multiple copies of the original 
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graph, each of which supports the hopping of one individ- 
ual boson species, and one obtains independent differential 
inequalities for the particle densities aj, s {t) = tr[hj tS p(t)] 
for each species. 

Beyond general bosonic models, our arguments also apply 
to fermions and mixtures of bosons and fermions [ 1 1 1, and in 
fact even to anyonic systems. Again, in a first step one can 
decouple the individual species of particles (which mutually 
commute) to hop on independent graphs. Then, it is easy to 
check that our arguments work independently of the statistics 
of the particles, since [fij , H] in Eq. evaluates to the same 
expression in terms of the fermionic (anyonic) creation and 
annihilation operators. Even better, fermionic and anyonic 
systems yield stronger bounds for the higher moments, and 
thus for the scenario of general local observables: In Eq. Q, 
can be bounded by 1 instead of 7V P_1 , which yields a 

bound otj(t) < CNq e vt ~ l on the higher moments. Corre- 
sponding results also follow for spin systems, as these can be 
described as hardcore bosons. 

Our arguments work not only for unitary theories, but also 
for certain types of dissipative (Markovian) models, extend- 
ing lfT2ll to bosonic systems. For instance, in the practically 
relevant case of a bosonic system with particle losses, we have 
that 

p(t) = -i [iW] - ({b%p(t)} ~ 2b jP (t)bfj . 

3 

Therefore, 

6>j(t) = -itr([hj,HBn]p(t)) - Xtx(hjp(t)) , 

which shows that the contribution from the dissipative term to 
ctj is negative; thus, tighter differential inequalities and thus 
a lower speed of sound than in the Hamiltonian case can be 
obtained. 

To conclude, we have proven that there is a speed limit for 
the propagation of information in a system of interacting par- 
ticles. This result is particularly relevant for the case of bosons 
on a lattice, as bosonic systems cannot be assessed using the 
established techniques of Lieb-Robinson bounds due to the 
unboundedness of the bosonic hopping operator. Our argu- 
ment applies equally to bosonic, fermionic, anyonic, and spin 
systems, as well as mixtures thereof, with arbitrary interaction 
terms between the particles, and can be generalized to also ad- 
dress systems with dissipation. 

The key point that allowed us to make statements about the 
propagation of information in bosonic systems beyond Lieb- 
Robinson bounds was first to focus on a subset of observables 
relevant to detecting the propagation of particles, namely the 
number of particles present at each site, and second to de- 
vise a closed system of inequalities bounding the evolution of 
their expectation values. This allowed us to reduce the prob- 
lem of characterizing the full dynamics of the system, which 
takes place in a superexponentially large Fock space, to sim- 
ply keeping track of the dynamics of a relatively small num- 
ber of parameters. This considerably reduced the complexity 



of the problem and gave rise to an exactly solvable worst-case 
bound. 

The idea of studying information propagation by restricting 
to a specific set of observables and investigating the result- 
ing worst-case differential equation can also be applied to the 
study of continous systems. This can be done either by tak- 
ing an appropriate continuum limit of a lattice model, or by 
directly considering a corresponding differential equation for 
the particle density which is continuous in space. 
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